Conformation dependence of molecular conductance: chemistry versus geometry 
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Recent experiments by Venkatamaran et al. [Nature (London) 442, 904 (2006)] on a series of molecular 
wires with varying chemical compositions, revealed a linear dependence of the conductance on cos^O, where 
9 is the angle of twist between neighboring aromatic rings. To investigate whether or not this dependence 
has a more general applicability, we present a first principles theoretical study of the transport properties of 
this family of molecules as a function of the chemical composition, conformation and the contact atom and 
geometry. If the Fermi energy E-p lies within the HOMO-LUMO gap, then we reproduce the above experimental 
results. More generally, however, if E-p is located within either the LUMO or HOMO states, the presence 
of resonances destroys the linear dependence of the conductance on cos^O and gives rise to non-monotonic 
behaviour associated with the level structure of the different molecules. Our results suggest that the above 
experiments provide a novel method for extracting spectroscopic information about molecules contacted to 
electrodes. 

PACS numbers: 85.65.+h,73.63.-b,73.40.Cg,73.40.-c 
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Single-molecule electronics poses many fundamental chal- 
lenges for chemistry, physics and engineering, partly because 
a molecule attached to metallic electrodes is a quantum ob- 
ject that lives between the traditional disciplines of chemistry 
and physics. Challenges also arise, because many electrical 
properties of single-molecule devices are exponentially sen- 
sitive to changes in the environment and to details of the 
contact to electrodes. As a consequence, experiments using 
mechanically-controlled break junction techniques (MCBJ) 
III] and scanning tunneling microscopy (STM) |2|, |3|] yield 
broad distributions of measured transport properties, corre- 
sponding to slightly different attachments of niolecules and 
in some cases, to environmental fluctuations Given 
these sensitivities, it is remarkable that systematic experi- 
mental studies can yield clear trends associated with generic 
molecular features. One such example is a study of the ef- 
fect of tilting the angle of contact between the long axis of a 
molecular wire and the plane of the locally-flat electrode 0]. 
Another example is a recent study of single-molecule electron 
transport as a function of molecular conformation 
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While both of the above studies reveal that geometry plays 
an important role in controlling electron transport through sin- 
gle molecules, the experiments of Ref. 7 contain an even 
stronger suggestion, namely that for the set of molecules stud- 
ied, changes in chemical composition are irrelevant. In these 
experiments, each molecule in the chosen series, (see Fig. 
([T]i), possessed two phenyl rings linked by a single carbon- 
carbon bond. Different chemical groups on the molecules 
caused the static angle of twist 6 to vary from molecule to 
molecule. The highest conjugation and therefore the highest 
conductance is expected when the two rings lie in the same 
plane (ie 6 = 0). However, when the rings are rotated relative 
to each other, the overlap between the tt orbitals decreases and 
the electrical conductance follows a cos^6' law. This simple 
dependence on geometry is remarkable, since different chemi- 
cal side groups were used to produce different angles. In other 
words, it appears that chemistry is irrelevant and geometry is 



FIG. I : GGA-relaxed configurations of all molecules studied in this 
work capped with NH2 . Grey, green, blue, red and orange vertexes 
correspond to H, C, N, F and CI atoms, respectively. 



everything. 

In this article, we aim to understand whether or not the 
dominance of geometry over chemistry can be widely ex- 
pected or if this is a peculiar property of the set of molecules 
measured. Close inspection of the measured conductances in 
Ref. 101 reveals that they do not follow a perfect cos^6' law 
and it is of interest to understand these deviations. The in- 
fluence of the chemical composition can be separated in two 
parts. On the one hand, different atoms or molecular groups 
can be attached to the aromatic rings to produce a desired an- 
gle. These include hydrogen atoms, alkanes, other aromatic 
chains, nitrogen, fluorine and chlorine atoms. Such groups 
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s calculated using LDA and GGA. 
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define the ground-state angle 0, which can vary between 
and almost 90 degrees. On the other hand, it is also possi- 
ble to employ different end atoms, which attach the molecules 
to the leads. In the experiments, Venkatamaran et al. Hi 
used molecules capped with either thiol (SH) or amino groups 
(NH2) and found that nitrogen termination gave more repro- 
ducible results. 

In what follows, we present a theoretical analysis of the 
structural configurations and transport properties of a series of 
molecules shown in Fig. ([T]i, with thiol or amino end-groups. 
Results were obtained using the SIESTA implementation of 
DFT |9], which employs norm-conserving pseudopotentials 
and linear combinations of atomic orbitals [10]. To treat ex- 
change and coiTelation we used the generalized gradient ap- 
proximation (GGA) [11], but we also made some tests with 
the local density approximation (LDA) [12], which allowed 
us to check the influence of the energy functional on struc- 
tural and transport properties ITsll . Transport properties were 
calculated with the SMEAGOL code (Q, which is interfaced 
to SIESTA and self-consistently computes the charge density, 
the electron-transmission coefficients and the I-V characteris- 



tics lUJj. 

To find the most stable structure we started from a twisted 
configuration, with an angle between and 90 degrees, and 
then allowed the molecule to freely relax. Since the pres- 
ence of local minima could have prevented the rings from ro- 
tating to the most stable configuration, we used various ini- 
tial angles to ensure that we reached the absolute minimum. 
9 = corresponds to a meta-stable configuration for most 
of the molecules. For other starting angles we found that the 
molecules always relaxed to the most stable configuration. Ta- 
ble dUl shows the ground-state angles calculated with LDA and 
GGA. Both functionals give similar results, although in most 
cases the LDA underestimates the angle. Both sets of results 
compare well with previous GGA values, computed in Ref. 
0. Table dUialso shows the dependence of 6 on the end-groups 
and reveals that the end group can change by up to 15%. 

We used the relaxed conformation of the pristine molecule 
as a starting point to compute the transport characteristics be- 
tween fee (111) gold electrodes. We assumed that the pro- 
tecting hydrogens of the amino and thiol groups were lost 
when the molecules were attached and the angle between the 
molecule and the surface was 90 degrees. We carefully deter- 
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FIG. 2: Transmission curves for molecules 1-8 between gold 
contacted by a surface gold adatom. 



leads 
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FIG. 3: Transmission curves for molecules 1-8 between gold leads 
contacted by a three-atoms hollow site. 



mined the distance between the molecule and the gold leads 
by studying the energy as a function of the distance of a sim- 
plified molecule composed of a benzene attached to only one 
N or S on top of a (111) gold slab made of 5 layers. With this 
distance and the molecule in its ground-state conformation, 
we calculated the transmission coefficient T{E) for electrons 
of energy E passing through the device. 

The transmission coefficients of all molecules between gold 
(1 1 1) electrodes computed with GGA are shown in Figs. (|2|l 
and (O 111 61] . We found that nitrogen termination hole-dopes 
the system, compared to the thiol termination, since the for- 
mer moves the HOMO closer to the computed Fermi energy 
Ep = 0. We also studied two possible contact geometries, 
namely a hollow position (Fig. ([3]l), where the contact atom 
sat on top of a 3-atoms hollow of the (111) surface, and a 
top position (Fig. (|2]i), where a gold adatom was placed on 
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top of the hollow position and in contact with the molecule. 
By comparing both contact configurations we found that the 
top position produces a much weaker coupling 11711 . mani- 
fested by the fact that the transmission coefficients in Fig. (|2]l 
exhibit a series of sharp Breit-Wigner resonances of height 
equal to unity. In the hollow position, Fig. (O shows that the 
Breit-Wigner description breaks down, because the coupling 
is much stronger and consequently, the resonances are merged 
and smeared out. The use of LDA and/or DZP changes the 
height and shape of the LUMO peak, but gives a similar value 
for the T{E) at the Fermi energy and almost the same HOMO 
peak. 

Close inspection of the structure of T{E) reveals informa- 
tion about the molecular states. Fig. Q shows the evolution 
of the peaks in the top configuration as 9 increases and re- 
veals that both HOMO peaks, which initially are clearly sepa- 
rated, move together and eventually merge for angles close to 
90 degrees. This behavior is reminiscent of bonding and anti- 
bonding states with a coupling matrix element, which depends 
on the overlap between the tt conjugated states of both rings. 
From the same figure, it is also clear that the width of the res- 
onances does not decrease as the angle changes, which means 
that the couplings to the leads (the F matrices) are the same for 
all molcules and the observed changes have an intra-molecular 
origin. In the hollow configuration. Fig. ^ shows that T{E) 
is characterized by the presence of two broad peaks associ- 
ated to HOMO and LUMO levels. As 6 increases, the width 
of the HOMO decreases, which is again due to the fact that 
the splitting between bonding and anti-bonding states, within 
the HOMO peak, decreases with increasing 6. 

Having examined T{E) for a range of energies, we now 
turn to the low-voltage, low-temperature conductance G{9) ~ 
GoT{Ef, 9), where Go = 2e^//i, Ep is the Fermi energy and 
for ease of notation we have now made the dependence of the 
transmission coefficient on 6 explicit. In figures (|2]l and (O, 
the zero of energy is chosen to coincide with the computed 
Fermi energy Ep. In an experiment, Ep may differ from Ep 
for a number of reasons, including the presence of a dielec- 
tric environment, such as air or water [0, Etl- Furthermore, 
even in the absence of environmental effects, differences can 
arise from the absence of self-interaction corrections in the 



exchange-correlation functional OlSll . For these reasons, we 
shall treat Ep as a free parameter and examine the 6 depen- 
dence of G for a range of Ep. 

Our first observation is that when the Fermi energy is inside 
the HOMO-LUMO gap, the conductance is approximately a 
linear function of cos^6', in agreement with the experiments of 
Ref. 0- This result is illustrated in Fig. ^ (a), which shows a 
comparison between experimental results for the normalised 
conductance G = G{0) / G(0) and theoretical results 1 19], ob- 
tained in the hollow configuration, for both N and S couplings. 
(In the experimental curve we also use our GGA-computed 
angular values.) As can be seen the, agreement with the ex- 
perimental curve is rather good. When LDA is used, slight 
differences in the computed the angle can shift the position of 
the points in the graph, but since the conductance values are 




FIG. 4: (a) Zero-bias conductance in the hollow configuration, for 
sulfur contact (circles), nitrogen contact (squares) and values from 
Ref. (triangles). All cases have been normalized to the 6 — 
value, (b) and (c) computed low-bias conductance G as a function of 
cos^O and the Fermi level for the the top configuration and nitrogen 
contact. The inset shows the curves of graph (b) normalized to the 
^ = value. 



also different at those angles, the dependence is still linear and 
the curve sits on top of those shown in Fig. (H^). 

Fig. (in (b) shows plots of T{Ep,6) obtained in the top 
configuration, for several values of Ep. Again we find a lin- 
ear dependence on cos^6', when the Fermi energy sits in the 
HOMO-LUMO gap. If these curves are normalized to the 
6* = value, as shown in the inset of Fig (HJ^), all of them 
sit on top of each other, which indicates the robustness of the 
linear dependence inside the HOMO-LUMO gap. 

In contrast with the above behaviour, when Ep approaches 
either the HOMO or the LUMO orbitals, the 9 dependence 
of the molecular resonances destroys the linear behavior and 
T{Ep, 9) can even become non-monotonic. This is illustrated 
in Fig. (IDl (c) for the case of N bonding to a top site. Similar 
behavior is found in the hollow configuration or when N is 
substituted by S as bonding atom. This implies that by gating 
1 2M the molecules, the measured 9 dependence could be tuned 
to be either linear or non-monotonic, depending on the value 

of £;f. 
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In summary, we have studied the transport properties of a 
series of molecular wires as a function of their conformal, 
chemical and contact configuration. We find that geometry 
is the dominant factor when the Fermi level sits inside the 
HOMO-LUMO gap, since different side groups do not af- 
fect the linear dependence. However, when the Fermi level 
approaches either the HOMO or LUMO resonances, chem- 
istry comes into play, because the chemical composition of 
the pristine molecules causes the positions of the HOMO and 
LUMO resonances to differ from molecule to molecule. This 
result shows that experiments such as those in Ref. tI yield 
important spectroscopic information about molecules in con- 
tact with electrodes, since the experimentally-measured linear 
dependence on cos^9, allows one to conclude that the Fermi 
energy in these experiments lies in the HOMO-LUMO gap. 
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